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SUMMARY 

Work completed during contract: 

1. Completed writing and editing of “The Middle and Upper Atmosphere of Neptune”, by J. 

Bishop, S. K. Atreya, P. N. Romani, G. S. Orton, B. R. Sandel, and R. V. Yelle, to appear 
as a chapter in the book Neptune and Triton^ D. Cruikshank and M. S. Matthews, eds.. 
University of Arizona Press. (In press, see Attachment 1). 

2. Completed work on UVS occultation lightcurve modeling using photochemical models fitting 

ground-based and Voyager IRIS infrared data. This work is documented in the attached 
DPS poster materials (Bishop et aJ. 1994) (Attachment 2). 

3. Supported the Voyager IRIS data modeling by Romani et al. (1993) (see Attachment 3). 
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TASKS ACCOMPLISHED 


While it was a distinct honor to have been asked to take the lead in compiling a review 
chapter, it was a very time-consuming task. I believe that we managed to structure the 
presentation so as to be suitable for persons entering this field of study {e.g., students and 
post-docs) while covering all the pertinent topics in sufficient detail to be useful as a technical 
review. The galleys were checked and corrected for the final time back in February 1995 (see 
Attachment 1). . 

Work on completing our analysis of the Voyager UVS solar occultation data acquired 
during Neptune encounter is essentially complete, as testified by the attached poster materials 
(Attachment 2). These materials provide a self-contained summary of our achievements in 
this area. A paper has been in preparation for some time, and requires only about a week’s 
- worth of effort- to complete; this will be returned to in December of this year, and submission 
of the paper to Journal of Geophysical Research - Planets should occur in January. The 
photochemical modeling addresses the recent revision in branching ratios for radical production 
in the photolysis of methane at H Lyman a implied by the lab measurements of Mordaunt 
et al. (1993). It should be noted that Mordaunt et al. did not present definitive branching 
ratios but instead simply demonstrated that methyl production occurs about half the time; the 
branching ratios for production of other radicals are still uncertain, so that it is necessary to 
ascertain the variation in modeling results associated with the Mordaunt et al. uncertainties. 
As it turns out, the modification of CH 4 branching ratios suggested by the Mordaunt et al. 
results (discussed in the poster materials) does not lead to significantly different photochemical 
modeling results compared with modeling incorporating the branching ratios advocated by 
Slanger and Black (1982). Yelle et al. (1993) carried out an independent analysis of the 
UVS occultation lightcurves using “objective” inversion techniques similar in spirit to the 
techniques we applied in our analysis of the Uranian solar occultation data (Bishop et al. 
1990). We had initially tried to apply these techniques to the Neptune data but found the 
data to be too noisy for reliable results, contrary to the claims of Yelle et al. Instead, we rely 
on the forward modeling approach described in Bishop et al. (1992), augmented with a more 
thorough incorporation of the Romani et al. (1993) photochemical model. While claiming to 
address photochemical issues, Yelle et al. did not carry out photochemical modeling calculations 
but rather relied on back-of-the-envelope estimates for column production and loss rates and 
product species fluxes (principally ethane, C 2 H 6 ) that axe highly questionable. The results 
of their analysis were forcefully presented; nevertheless, it has been necessary to assess their 
conclusions carefully to avoid placing too much reliance upon them. 

The software generated in this effort hcis been useful for checking the degree to which 
photochemical models addressing other datasets (mainly infrared) are consistent with the 
UVS data. This is documented in the Romani et al. (1993) paper -(A-btachment 3). This 
work complements the UVS modeling results in that the IR data refer to deeper pressure 
levels {p > 0.1 mbar); as regards the modeling of UVS data, the most significant result is 
the convincing support for the presence of a stagnant lower stratosphere. Evidence for strong 
dynamical (mixing) transport of minor constituents at shallower pressures is provided by the 
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UVS data analysis, discussed in the Attachment 2 materials. 


All aspects of this work were carried out in close collaboration with Dr. Paul N. Romani 

(NASA/GSFC). 
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Observations peitaining to the middle and upper atmosphere of Neptune obtained &om 
Earth have, until lecently, been limited to hydrocarbon thermal emission spectra and 
refiractive stellar occultations. The Voyager 2 encounter provided our Srst detailed 
look at these atmospheric regions. RSS and IRIS measurements have confirmed the 
p — T structure near the tiopopause inferred from pre-Voyager infrared measurements 
but have also revealed a steeper temperature gradient at smaller pressures and pro- 
nounced latitudinal variations, with tropopause temperatures in the range 50 K to 
57 K. The pre-Voyager estimate of a stratospheric methane mixing ratio of 2%. more 
than 2 otdm of magnitude in excess of the tropopause cold trap value, has been 
revised downward considerably (10~*-I0~^) but oversatuiation is still indicated. In 
contrast. Voyager measurements of stratospheric CjHj and CjHt abundances are in 
very close agreement with estimates from Earth-based observations: x 10~* aird 

X 10”*. respectively, near 03 mbar (T w 160 K). Vertical transport in the upper 
stratosphere is fairly vigorous, with eddy mixing coefficient values in the range 3 to 
10x10* cm' s~' near the 02 ^bar level as determined from the U VS solar occultation 
data. Eddy mixing in the lower stratosphere {pZ2 mbar) is sluggish by comparison, 
with values on the order of 1 to 3 x 10^cm^~' indicated by photochemical modeling 
of the infiared emission spectra. Submiilimeter emissions by HCN and CO in the 
suatosphere, first detected in 1991. have raised new questions regarding stratospheric 
chemistry and transport on Neptune; modeling of the abundances of these species, in 
conjunction with continued photochemical modeling of hydrocarbon thermal emis- 
sions and the UVS solar occultation data, is likely to lead to a characterization of 
eddy mixing throughout the stratosphere and to tighter constraints on stratospheric 
N: and CHt mixing ratios. Other prominent modeling issues of current interest are 





m 


ATTACHMENT 1 



ATTACHMENT 2 


PHOTOCHEMICAL MODELING OF VOYAGER UVS SOLAR 
OCCULTATION LIGHTCURVES ACQUIRED AT NEPTUNE 


James Bishop (Computational Physics Inc.) 

Paul N. Romani (NASA/Goddard Space Flight Center) 

Sushil K. Atreya (University of Michigan) 

Ingress and egress occultations of the Sun by Neptune were recorded with 
the UVS instrument on Voyager 2 [Broadfoot et al. 1989]. Previous model- 
ing of the 125-140 nm lightcurves indicated eddy mixing coefficient values of 
3 — 10 X 10® cm^ s“^ near the 0.2 /xbar level and methane mixing ratios in the 
lower stratosphere on the order of 1 — 3 x lO""^ [Bishop et al. 1992]. These re- 
sults should be insensitive to photochemical details, provided methane is the 
main source of opacity at these wavelengths. The UVS lightcurves at longer 
wavelengths (140-153 nm) are expected to be dominated by the opacity of 
C 2 species (ethane, acetylene, ethylene) and perhaps higher order organics. 
At still longer wavelengths, H 2 Rayleigh scattering is a major opacity source. 
We will present photochemical models giving good agreement with both sets 
of UVS lightcurves at wavelengths 125-165 nm and with C 2 H 6 and C 2 H 2 
abundances near 0.5 mbar derived from IRIS mecisurements [Bezard et al. 
1991]. The current photochemical model incorporates several updates, in- 
cluding the recent revision in CH 4 photolysis branching ratios at Lyman a 
[Mordaunt et al. 1993]. These model fits strongly suggest a region of en- 
hanced eddy mixing near the 10 /ibar level (K ~ 10® cm^ s~^), decreasing at 
higher altitudes. Comparisons between egress (~ 49°S latitude) and ingress 
(61°N) UVS results will be discussed, as well as the impact of uncertainties in 
mean stratospheric temperature and key reaction rates. In line with our ear- 
lier work, methane mixing ratios on the order of 10 “'^ are required to obtain 
good agreement between the photochemical models and the UVS lightcurves. 


This work supported by NASA Contract NASW-4806. 



UVS LIGHTCURVES 


The UVS solar occultation lightcurves relevant to stratospheric investigations fall into three 
groups: for Ac ~ 138 nm, the lightcurves exhibit half-light points in the 0. 1-0.3 //bar region, 
those with 138 ~ Ac ~ 153 nm probe progressively deeper into the stratosphere, while the 
half-light points of Ac ~ 153 nm lightcurves are at pressures ~ 0.1 mbar. In the figure below, 
a subset of the UVS solar occultation data used in our modeling is shown as transmission 
lightcurves K)/h{K), where z± is the minimum tangent ray height) averaged over 

0.96 s intervals and ordered according to the channel center wavelengths Ac at 500 km MTRH. 
Altitude resolutions are 9.5 km (ingress) and 12.3 km (egress), and the spectral resolution is 
~2.5 nm. The apparent solar diameter for both occupations (~4 km at ingress, ~20 km at 
egress) is small relative to the MTRH range over which total opacities drop from 0.9 to 0.1 
and can be ignored. Other aspects of the occultation (e.g., channel center offsets as functions 
of MTRH stemming from the limit cycle motion, convolution of model transmission spectra 
at 0.1 nm resolution with the UVS slit function) are incorporated in our modeling using the 
information provided to us by B. Sandel and R. Vervack (Univ. of Arizona). 


In addition to measurement noise, the data are afflicted by two main sources of error: internal 
instrument scattering, primarily of the solar Lyman a line, and pointing errors which introduce 
uncertainties in the ratioing. These errors and the derivation of a corresponding standard 
deviation ctd are discussed by Yelle et al. [1993]. In recognition of this, we display the UVS 
data as data ranges defined by /(^x, Ac)//o(Ac) ± (Td when making comparisons with model 


lightcurves. 
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IRIS Results: The IRIS measurements of stratospheric hydrocarbon emissions analyzed by 
Bezard et al. [1991] yielded acetylene and ethane mixing ratios of ~4 x 10“® (0.2 mbar) and 
~1.3 X 10~® (0.7 mbar), respectively; these results are based on a large selection of low spatial 
resolution data with field-of-view centers ranging from 10°S to 50°S latitude. Since the UVS 
egress occultation latitude lies within this range, the IRIS-derived mixing ratios have been 
used to help identify successful photochemical models under egress conditions. 




MODELING PROCEDURE 


The model lightcurves are constructed using the H2 density distributions provided by the 
model atmospheres and hydrocarbon (CH4, C2H2, C2H4, C2H6) density distributions obtained 
from solving the 1 -D continuity equation 


d$. 

dz 


= Pi - Li 


for each species i assuming steady-state conditions. (Variables are defined below). In terms of 
the mixing ratio for species i, the flux is 

= -JV(A + - DiNMHr' - 

Once the hydrocarbon abundances are evaluated, simulated UVS lightcurves are given by 


/*Ac+A«A foo V* 

^c)=j R{>^ - K) ■ /o(A) • exp -2 f dr • — — ■ n{z) ■ ^ fi{z)( 7 i{\) 

JXc — Aa^ L V — 'X/ ' i 


The main fitting parameters are the mixing ratio of methane in lower stratosphere (/t(CH 4)) 
and the parameters defining the eddy mixing coefficient profile. We have explored several 
different types of K profiles. In seeking out successful model fits, our first step is to identify 
(/t(CH4), Kii2{\c < 140 nm)) pairs providing the best agreement with the UVS data at 
wavelengths shortwards of 140 nm; K\j2{^c) is the eddy mixing coefficient at the half-light point 
of channel Ac- As previously mentioned, these lightcurves are expected to be least sensitive to 
photochemical modeling details. Using the constraints on /t(CH 4) and the strength of eddy 
mixing at the top of the region probed by the UVS lightcurves, the K profile is then varied (as 
allowed by the adopted form) to try to obtain agreement at longer wavelengths, which probe 
progressively deeper pressures. 


Variable Definitions 
fi 

Pi 

Li 

Di 

K 

Hi 

Hatm 


Z, r 


zx, rx 

A 

Ac 


A,A 

R{\ - Xc) 
lo 


<^i 

n 


mixing ratio of species i 

vertical flux of species i (molecules cm“^ s~^ ) 

chemical production rate of species i (molecules cm“^ s~^) 

chemical loss rate of species i (molecules cm“^ s~^) 

molecular diffusion coefficient of species i for H2~He atmosphere (cm^ s~^) 

atmospheric eddy mixing coefficient (cm^ s“^) 

partial pressure scale height of species i 

pressure scale height of the background atmosphere 

altitude, radial distance from center of planet 

minimum tangent ray height for line of sight to center of solar disk 
wavelength 

channel center wavelength at Zx 
UVS channel spectral width 

normalized instrument response at A in channel with center wavelength A^ 
solar flux at top of atmosphere (photons cm~^ s“^ A“^) 
photoabsorption cross section of species i (cm^) 
atmospheric number density at altitude z (cm“^) 



ATMOSPHERE MODELS 


In our modeling, we adopt a mean molecular weight of 2.393 [Conrath et al. 1991] with a helium 
mixing ratio /(He) of U.19. .Atmospheric pressure-temperature (p — T) profiles at pressures 
deeper than ~ 2 mbar at the ingress and mean egress latitudes are taken from Lindal [1992]. 
These are extended to 1.0 ^bar using a profile of the form presented by Hubbard et al. [1987]; 
after scaling to /(He) = 0.19, their analysis of the 20 August 1985 stellar occultation central 
flash indicates a temperature of 143 K near 0.42 mbar. The nominal 1 ^bar temperature is 
168 K, taken from Orton et al. [1992]. The high altitude segment is based on the preliminary 
analysis of the UVS ingress solar occultation data presented by Broadfoot et al. [1989]. 

Since the p~T profile is not well-constrained by Voyager flyby data at pressures in the range 
2 mbar to 0.01 /<bar, “warm” and “cool” models have been used in the modeling. These are 
defined by shifting the 1 //bar temperature of Orton et al. by ±20 K as suggested by numerous 
stellar occultations [Roques et al. 1994]. In these models, the RSS profiles and Hubbard et al. 
0.42 mbar point are not modified. 

The 1 bar radii at the UVS occultation latitudes, which are needed to place the p — T profiles 
on an altitude scale, are taken from Lindal [1992]: 24443 km at the ingress latitude of 61°N and 
24535 km at the mean egress latitude of 49°S. The altitude scale on the right hand ordinate 
axis in the figure below corresponds to the nominal ingress p — T model. 

The region where photochemically produced hydrocarbons are predicted to condense is also 
indicated [Moses et al. 1992; Romani et al. 1993]. 
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OTHER INPUTS 


Solar Irradiance Spectrum 

Photolysis rates are calculated using solar irradiance fluxes summed over 10 A intervals span- 
ning the wavelength range 1000-2500 A. At wavelengths longward of 1204 A, the fluxes are 
from the 1985 Spacelab-2 flight of SUSIM [see Bruckner et al. 1993], and a scaled (Fio.? = 207) 
Hinterreger spectrum [Hinteregger et al. 1981] is used at shorter wavelengths. The composite 
spectrum at 1 AU is shown below and is representative of solar maximum conditions; the 
line-integrated Lyman a flux is 6.5 x 10^^ ph cm~^ s“^. At Neptune’s distance from the Sun, 
the pervasive Lyman a skyglow associated with the local interstellar medium is not negligi- 
ble and has been included in the modeling; based on near-encounter LISM measurements by 
the Voyager UVS, the contribution to the CH 4 photolysis rate is 3 x 10“® s“^ (at the top of 
the atmosphere). The SUSIM spectrum at 1 A resolution is also employed in the lightcurve 
modeling. 

In our modeling, we adopt “local” conditions, i.e., we adopt model atmospheres and evaluate 
photolysis rates corresponding to conditions at the occultation latitudes at the time of en- 
counter. Our modeling of the UVS lightcurves suggest rapid eddy mixing at pressures 1 /ibar 
~ p ~ 0.1 mbar, so that the timescale for mixing (Tmix = H^/K, where Hn is the local atmo- 
spheric density scale height and K is the eddy mixing coefficient) is significantly smaller that a 
solar cycle (and consequently much shorter than a Neptune season). The alternate analysis of 
the UVS data by Yelle et al. [1993] suggests K w 10® cm^ s“^, which is smaller thaja the values 
of K we infer at the same pressures; even with this value, though, the mixing timescales are 
significantly smaller than seasonal or solar cycle durations. The effective solar zenith angles 
adopted in the ingress and egress modeling are 89.2° and 40.2°, respectively, corresponding to 
diurnally averaged solax illumination conditions at the time of encounter. 
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Photoabsorption Cross Sections 


The photoabsorption cross sections used in the opacity modeling axe shown in the accompa- 
nying figure, at roughly 2 A resolution. These are taken from: 

CH4 Lee and Chiang [1983] & Mount et al. [1977] 

C 2 H 2 Suto and Lee [1984] h Wu et al. [1989; personal communication] (at T = 155 K) 
C 2 H 4 R. Wu, personal communication (T = 155 K) 

C 2 H 6 Mount and Moos [1978] 

C 4 H 2 dicker and Okabe [J. Phys. Chem. 91:437, 1987] 

Superimposed are the cross sections convolved with the UVS instrument slit function. The 
channel center wavelengths for the ingress occultation data at a tangent altitude of 500 km 
axe also indicated. 

The photolysis channels and branching ratios included in the photochemical modeling are listed 
on the accompanying table. The recent revision in the direct quantum yield of CH3 from the 
photolysis of methane at Lyman a has been folded into our modeling. However, since the story 
of CH4 Ly-a photolysis is still incomplete, we have attempted to make reasonable estimates of 
the branching ratios for each energy-allowed channel. These estimates are discussed in detail 
on the accompanying page. 

Chemical Reaction Rates 

The chemical reactions included in the modeling are listed on the accompanying table along 
with the adopted rates. There have been several updates to the modeled chemical scheme 
since the study by Romani et al. [1993]. The most notable changes axe (1) the introduction of 
the photolysis channel CH4 + hi/ —* CH3 + H using the results of Mordaunt et al. [1993] and 
(2) revision of the chemical pathways involving the vinyl radical (C2H3) (see talk by Romani 
et al., this meeting). The chemical scheme is relatively complete for the C 2 hydrocarbons 
(see figure). Expanding the chemical scheme to a more complete coverage of reaction paths 
involving C3 and C4 species, similar to the photochemicaJ model presented by Moses et al. 
[1992], will be carried out in the near future. 

Condensation 

In the lower stratosphere, temperatures are cold enough to result in very large supersaturations 
(5) for several hydrocarbon species (principally C 2 H 6 and C 2 H 2 ) if these remained in the gas 
phcise. Consequently, a condensation loss is included. Our handling of this loss mechanism is 
discussed in detail in Romani et al. [1993]; the nominal placement for removal from the gas 
phase is dictated by the condition S > 1. The study of the relative efficiency of condensation 
processes by Moses et al. [1992] suggests that the condensation levels of product hydrocarbons 
should lie deeper in the atmosphere than indicated by the 5 > 1 condition; however, both 
the IRIS and UVS hydrocarbon data refer to pressure levels {p w 0.5 mbar and p ~ 0.1 
mbar, respectively) well above the condensation levels so that modeling of these data are only 
indirectly affected by the exact placement of these levels. The important point is the existence 
of an efficient loss mechanism acting at a deeper level that establishes downward fluxes of 
acetylene and ethane in the pressure regions sounded by the IRIS and UVS data. 






PHOTOLYSIS CHANNELS AND REFERENCES 



REACTION 


BRANCHING RATIOS 

REFERENCES 

Jla 

CH4 + hi/ 

+ H + H 

9(121.6) = 0.212 

*** see text *** 

Lee & Chiang (1983) J. Chem. Phys. 78:688 
Mount et al. (1977) Astrophys. J. 214;L47 

Jib 

CH4 + hj/ 

'CH2 -b H2 

9(121.6) = 0.282 

McNesby & Okabe (1964) Adv. in Photochem. 3:157 
Calvert and Pitts (1966) Photochemistry, J Wiley & Sons, Inc. 

Jlc 

CH4 + hj/ ^ 

CH + H + H2 

9(121.6) = 0.100 

Rebbert & Ausloos (1972) J. Photochem. 1:171 
Gorden & Ausloos (1967) J. Chem. Phys. 46:4823 

Jld 

CH4+hi/ ^ 

CH3 + H 

9(121.6) = 0.406 

Slanger & Black (1982) J, Chem. Phys. 77:2432 
Mordaunt et al. (1993) J. Chem. Phys. 98:2054 

Jle 

CH4 + hi/ 

ICH2 + H + H 

9(121.6) = 0.0 

Laufer k McNesby (1968) J. Chem. Phys. 49:2272 

J2a 

C2H2 "b h.u — > 

C2H + H 

9(A < 153) = 0.30 
9(153 < A < 186) = 0.12 
9(186 < A < 201) = 0.21 

Okabe (1981) J. Chem. Phys. 75:2772 
Okabe (1981) J. Chem. Phys. 78:1312 
Shin & Michael (1991) J. Phys. Chem. 95:5864 
McDonald et al. (1978) Chem. Phys. 33:161 

J2b 

C2H2 "b ht' — > 

C2 *b H2 

9(A < 201) = 0.10 

Suto & Lee (1984) J. Chem. Phys. 80:4824 
Wu et al. (1989) J. Chem. Phys. 91:272 
Nakayama k Watanabe (1964) J. Chem. Phys. 40:558 
R. Wu (1991) personal communication 

J3a 

C2H4 + hi/ 

C2H2 + H2 

9(A < 190) = 0.51 

Xia et al. (1991) J. Quant. Spect. Rad. Transfer 45:77 
R. Wu (1991) personal communication 

J3b 

C2H4 + hi/ 

C2H2 + H + H 

9(A < 190) = 0.49 

Kara k Tanaka (1973) Bull. Chem. Soc. Japan 46:3012 
Sauer k Dorfman (1961) J. Chem. Phys. 35:497 
Back k Griffiths (1967) J. Chem. Phys. 46:4839 






C2H6 

+ 

hi/ 

C2H2 + H2 4 - H2 

9(121.6) = 

0.25 





9(A < 170) 

= 0.27 

C2H6 

+ 

hu ^ 

C2H4 + H + H 

9(121.6) = 

0.30 





9(A < 170) 

= 0.14 

C2H6 

+ 

hi/ — »• 

CH4+ 'CH2 

9(121.6) = 

0.25 





9(A < 170) 

= 0.02 


9(121.6) = 0.13 
q{\ < 170) = 0.56 

9(121.6) = 0.08 
q{\ < 170) = 0.01 

q{\ < 166) = 0.20 

9(A < 166) = 0.10 
9(166 < A < 207) = 0.06 

J5c C4H2 + h// — > C2H -f- C2H 9(A < 166) = 0.03 

9(166 < A < 207) = 0.01 

J5d C4H2 + hi/ ^ C4H;* 9(A < 166) = 0.67 

9(166 < A < 207) = 0.93 
9(A > 207) = 1.00 


J 4 d 

C2He + hi/ ^ C2H4 + H2 



J 4 e 

C2H6 + hi/ — 

CH3-I-CH; 

J 5 a 

C4H2 + hi/ - 

C4H + H 

J 5 b 

C4H2 + hi/ - 

-> C2H2 + C2 


Mount & Moos (1978) Astrophys. J. 224:L35 

Calvert & Pitts (1966) Photochemistry, J. Wiley & Sons, Inc. 

Akimoto et al. (1965) J. Chem. Phys. 42:3864 
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Hampson & McNesby (1965) J. Chem. Phys. 42:2200 
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REACTIONS, RATES AND REFERENCES 



REACTION 

RATE EXPRESSION 

REFERENCES 

R1 

'CH 2 + H 2 ^ CH 3 + H 

k = 9.24 X 10'^^ 

Allen et al. (1992) Icarus 100:527 
Langford et al. (1983) J. Chem. Phys. 78:6650 
Ashfold et al. (1981) Chem. Phys. 55:245 
Baun et al. (1970) J. Chem. Phys. 52:5131 

R2 

'CH 2 + CH 4 -> CH 3 + CH 3 

ik = 6.0 X 10 -^^ 

Bohland et al. (1985) Ber. Bunsenges. Phys. Chem. 89:1110 

R3 

CH + CH 4 ^ C 2 H 4 + H 

A: = 5.0 X 10-'^ 
kmax = 1.7 X 10-'° 

Baulch et al. (1992) J. Phys. Chem. Ref. Data 21:411 
upper limit is T = 167 K rate from Berman & Lin (1983) 
Chem. Phys. 82:435 

R4 

CH + H 2 ^ CH 3 

Ajo = 8.75 X 10-°' 
k^ = 8.3 X 10-" 

h — /CooAToA/ 

~ kco+koM 

fit to data in Berman & Lin (1984) J. Phys. Chem. 81:5743 

R5 

CH 3 + H CH 4 

ko = 4.0 X 10-°° 
k^ = 4.7 X 10-'° 

L _ kc^kpM , p 
~ koo+koM ^ 

with fc = 0.902 - 1.03 X 10-° 

Brouard et al. (1989) J. Phys. Chem. 93:4047 
T and /„ = 1 

R 6 

CH 3 + CH 3 ^ C 2 H 6 

ifc<, = 3.5xl0-"r-"°°e-'°°°/^ 
k^ = 6.0 X 10-" 

Baulch et al. (1992) J. Phys. Chem. Ref. Data 21:411 
N.B. bath gas is Ar 


~ koo+koM ^ 

with fc = 0.38e"^/^3 + 
and fn = 0.75 - 1.27 logjo/c 










R7 

H + H ^ 

H 2 

jfc = 2.7 X 10-=^^ T-° ® 

R8 

H + C 2 H 2 

^ C 2 H 3 

h = 3.3 X 10-^° 

ifcoo = 1.4 X 10-11 e-i3°°/^ 

L _ JkajksiM.. . p 
kco+kcM 

with fc = 0.44 and /„ = 0.75 - 

R9 

H + C 2 H 3 

— > C 2 H 2 + H 2 

fc = 5.0 X 10-11 

RIO 

H 2 + C 2 H 3 

C 2 H 4 + H 

k = 5.01 X 10-20 

Rll 

H + C 2 H 4 

^ C 2 H 5 

ko = 2.15 X 10-29 

k^ = 4.95 X 10-11 e-1051/^ 

h — kookoM 
~ kcc+koM 

R12 

H + C 2 H 5 

CH 3 + CH 3 

k = 7.95 X 10-11 6-12^2’ 

R13 

^CH 2 + H 

^ CH 3 

ifco = 3.1 X 10-^° 
koo = 1.5 X 10-1° 

— kookoM 

~ fco=+fcoM 

R14 

CH 3 + ^CH 2 ^ C 2 H 4 + H 

k = 7.0x 10-11 


R15 ^CH 2 + C 2 H 2 


Products 


^ = 1.99 X 10“^^ g-3332/T 


Baulch et al. (1992) J. Phys. Chem. Ref. Data 21:411 
N.B. bath gas is H2 

Baulch et al. (1992) J. Phys. Chem. Ref. Data 21:411 
N.B. bath gas is He 

1.27 logio/c 

see R64 

Tsang & Hampson (1986) J. Phys. Chem. Ref. Data 15:1087 
Fahr et al. (1994) in preparation 

Lightfoot & Pilling (1987) J. Phys. Chem. 93:3373 
(fast rate) 

Pratt & Wood (1984) J. Chem. Soc. Faraday Trans. 80:3419 
Gladstone (1983) Ph.D. Thesis, Cal. Inst, of Tech. 

Tsang & Hampson (1986) J. Phys. Chem. Ref. Data 15:1087 
Bohland et al. (1986) Symp. Int. Combust. Proc. 21:841 











R17 C 2 H + H 2 


C 2 H 2 + H 


)fc = 1.8x 10-“e-i°9o/r 


R18 

C 2 H + CH 4 -> C 2 H 2 + CH 3 

= 1.7x 10-^'e-5«/^ 

R19 

C 2 H + C 2 H 6 ^ C 2 H 2 + C 2 H 5 

IS: = 3.6 X 10"^^ 

R20 

C 2 H + C 2 H 2 ^ C 4 H 2 + H 

A: = 1.1 X 10-^°e28/^ 

R21 

P-C 3 H 4 + H ^ CH 3 + C 2 H 2 

k = 9.62 X 10‘^2 g-is 6 o/T 

R22 

'CH 2 + H 2 ^ 3 CH 2 + H 2 

k = 1.26 X 10"^^ 

R23A 

C 2 H 3 + C 2 H 3 ^ C 2 H 2 + C 2 H 4 

A: = 2.4 X 10-^^ 

R23B 

C 2 H 3 + C 2 H 3 ^ C 4 H 6 

ko = 1.3 X 10-22 
A;«, = 1.2 X 10-^“ 

h — knokoAf 

~ koo+koM 

R24 

C 4 H + H 2 C 4 H 2 + H 

A:24 = ki7 

R25 

C 4 H + CH 4 C 4 H 2 + CH 3 

k25 = A:i 8 

R26 

C 4 H + C 2 H 6 C 4 H 2 + C 2 H 5 

k26 = A:i9 

R27 

C 4 H + H ^ C 4 H 2 

k27 = ki6 

R28 

C 4 H + C 2 H 2 ^ C 6 H 2 + H 

k28 — ^20 


Tsang & Hampson (1986) J. Phys. Chem. Ref. Data 15:1087 


Koshi et al. (1992) J. Phys. Chem. 96:9839 

Opansky et al. (1993) 5th Inti. Cnf. on Lab. Research 
for Planetary Atmospheres, Boulder CO (poster Ch-01) 


Lander et al. (1990) J. Phys. Chem. 94:7759 


Pedersen et al, (1993) J. Phys. Chem. 97:6822 


Wagner & Zellnar (1972) Ber. Bunsenges. Phys. Chem, 76:518 


see R1 


Fahr et al. (1991) J. Phys. Chem. 95:3218 


ko'. Laufer et al. (1983) Icarus 56:560 

koo' Fahr et al. (1991) J. Phys. Chem. 95:3218 


in lieu of firm evidence, estimating that C 4 H rates are the 
same as their C 2 H analogues, see: 

Kiefer et al. (1990) Int. J. Chem. Kin. 22:747 
Tanzawa Gardiner (1980) J. Phys. Chem. 84:236 
Frank & Just (1990) Combust. Flame 











R29 

R30 

C 4 H + C 4 H 2 ^ C 8 H 2 + H 
C 2 H + C 4 H 2 -> C 6 H 2 + H 

^29 — ^20 
^30 = ^20 

R31 

C4Hr + C4H2 C 8 H 2 + H 2 

Jfc = 1.5 X 10-^° 

R32 

C 4 H 2 + H ^ C 4 H 3 

fco = 1.0 X 10 - 2 ® 

itoo = 1.39 X i0-^o^-ns4/T 

7„ — kookoM 
~ fcoo+fcoM 

R33 

C 4 H 3 + H — > C 4 H 2 + H 2 

A: = 1.2 X 10-“ 

R34 

C 4 H 3 + H -> C 2 H 2 + C 2 H 2 

fc = 3.3 X 10-“ 

R35 

C 2 + H 2 C 2 H-HH 

k = 1.77 X 10-10 

R36 

C 2 + CH 4 ^ C 2 H + CH 3 

k = 5.05 X 10““ 

R37 

C 2 H + C 2 H 4 -> C 4 H 4 + H 

A; = 1.3 X 10-1° 

R38 

C 2 H 5 + C 2 H 5 C 2 H 6 + C 2 H 4 

A: = 2.4 X 10-12 

R39 

3CH2 + ^CH 2 ^ C 2 H 2 + H + H 

A: = 3.1 X 10-1° 

R40 

1 CH 2 + CH 4 ^ 3 CH 2 + CH 4 

A: = 1.2 X 10-“ 


R41 CH 3 + C 2 H 3 CH 4 + C 2 H 2 


k = 3Ax 10"“ 


in lieu of firm evidence, estimating that C4H + C 4 H 2 
and C 2 H + C 4 H 2 have the same rate as R20 
(see references for R24-R28) 

estimating to be kinetic rate limit 

Nava et al. (1986) J. Geophys. Res. 91:4585 
N.B. ka is an estimate and also a lower limit 

Yung et al. (1984) Astrophys. J. Supp. 55:465 
Yung et al. (1984) Astrophys. J. Supp. 55:465 
Pitts et al. (1982) Chem. Phys. 68:417 
Pitts et al. (1982) Chem. Phys. 68:417 
Lander et al. (1990) J. Phys. Chem. 94:7759 
Baulch et al. (1992) J. Phys. Chem. Ref. Data 21:411 
Darwin (1989) Ph.D. Thesis, Univ. of Calif. Berkeley 
see R2 

Fahr et al. (1991) J. Phys. Chem. 95:3218 












R42 

CH 3 + C 2 H 5 — > CH 4 + C 2 H 4 

k = 1.28 X 10 -^^ T-°-32 

R43 

C 2 H 5 + H C 2 H 6 

ito = 5.5 X 10-23 

ifcoo = 1.5 X 10-^3 g-440/T 
7 . _ koakaM 
~ fcoo+fcoA<f 

R44 

CH 4 + C 2 H 3 CH 3 + C 2 H 4 

A: = 2.4x 10-24 T4 02e-2754/T 

R45 

^CH 2 + H CH + H 2 

ifc = 4.7 X 10-1° e-3™/^ 

R46 

C 2 H 3 + C 2 H 5 ^ I-C 4 H 8 

fc = 1.9 X 10-2^ 
A: = 2.5 X 10-11 

1 . _ knokoM 

^ ~ kcc+kcM 

R47 

C 2 H 3 + C 2 H 5 ^ C 2 H 4 + C 2 H 4 

k = 2.81 X 10-12 

R48 

C 2 H 3 + C 2 H 5 ^ C 2 H 2 + C 2 H 6 

II 

00 

R49 

3 CH 2 + C 2 H 5 CH 3 + C 2 H 4 

A: = 3.0 X 10-11 


R50 + C 2 H 3 


CH 3 + C 2 H 2 A: = 3.0 X 10 -" 


Tsang & Hampson (1986) J. Phys. Chem. Ref. Data 15:1087 
Baulch et al. (1992) J. Phys. Chem. Ref. Data 21:411 
estimating that ki 2 = 0.04 • Aqo from R60 

Gladstone (1983) Ph.D. Thesis, Cal. Inst, of Tech. 


Tsang & Hampson (1986) J. Phys. Chem. Ref. Data 15:1087 

Zabarnick et al. (1986) J. Chem. Phys. 85:4373 

ko\ Laufer et al. (1983) Icarus 56:560 

koo- Tsang h Hampson (1986) J. Phys. Chem. Ref. 

Data 15:1087 

estimating that the disproportion to recombination ratio is 
0.23 (average of C 2 H 3 and C 2 H 5 self reactions, 0.3 and 0.15 
respectively, from Tsang h Hampson (1986) J. Phys. Chem. 
Ref. Data 15:1087) and taking k^^ = k^^ 

Tsang h Hampson (1986) J. Phys. Chem. Ref. Data 15:1087 

Tsang h Hampson (1986) J. Phys. Chem. Ref. Data 15:1087 








R51 

C2H5 + C2H5 ^ C4Hio 

ko = 6.59 X 10-® T- 6.39 e- 301 /T 
koo = 1.9 X 10-^1 

h — kcakoM 
~ kocy+koM 

R52 

CH + C2H2 ^ C3H2 + H 

A: = 3.5 X 10-^° 

kmax — 5-3 X 10“^° 

R53 

CH + C2H4 ^P-C3H4 + H 

k = 2.2 X 10~^° 

R54 

CH + C2H4 ^ U-C 3 H 4 + H 

kmax = 7.1 X 10-'° 
F = 0.5 
k 5 z = F • k 
k 54 = {1 - F) • k 

R55 

CH + C2H6 —* Products 

k = 1.8 X 10-'° 
kmax = 4.4 X 10-'° 


R56 P-C 3 H 4 + H ^ C 3 H 5 ko = 1.636 x 

fcoc = 6.0 X 10-11 e-1233/r 

L feopfegAf 

“ k^+koM 


R57 C 2 H 3 + C 2 H 5 


CH3 + C3H5 k — 2.5 X 10 11 — A;46 


ko'. Laufer et al. (1983) Icarus 56:560, as reported in 
the NIST Chemical Kinetics Database V5.0. 

N.B. Laufer et al. presented tabular data; the NIST 
equation reproduces the Laufer et al. values to ~ 1% 
koo' Baulch et al. (1992) J. Phys. Chem. Ref. Data 21:411 

Baulch et al. (1992) J. Phys. Chem. Ref. Data 21:411 
upper limit is T = 171 K rate from Berman et al. (1982) 
Chem. Phys. 73:27 

Baulch et al. (1992) J. Phys. Chem. Ref. Data 21:411 
upper limit is T = 160 K rate from Berman et al. (1982) 
Chem. Phys. 73:27 

N.B. rate is split between the two isomers by factor F. 

Baulch et al. (1992) J. Phys. Chem. Ref. Data 21:411 
upper limit is T = 162 K rate from Berman & Lin (1983) 
Chem. Phys. 82:435 

Whytock et al. (1976) J. Chem. Phys. 65:191 

Wagner Sz Zellnar (1972) Ber. Bunsenges. Phys. Chem. 76:518 

Tsang &: Hampson (1986) J. Phys. Chem. Ref. Data 15:1087 
estimating that ^43 + ^57 = fcoo from R46 






R58 CH 3 + C 2 H 3 -> C 3 H 5 + H 
R59 CH 3 + C 2 H 3 ^ C 3 H 6 


R 60 CH3 + C2H5 C3H8 ko = 1.01 X 10^° e-^s 97 /T 

= 3.2 X 10 -^° r -°-32 

U feopfep A/ 

^ “ fcoo+fcoM 


ko = 1.3 X 10-22 

iboo = 1.2 X 10-^° 
^59 - 

ks& — 0.002 • kao 


R61 CH + H 2 -> ^CH 2 + H k = 2.38 x 10 “'° 

R62 H + C 2 H 5 C 2 H 4 + H 2 A: = 3.0 X 10"^2 

R63 CH + H C + H 2 A: = 1.4x10-11 

R64 H + C 2 H 3 C 2 H 4 ko = 1.49 X 10-2^ 

koo = 1.55 X 10-1° 

L — kookpAf 
““ kco+koM 


ko'. Laufer et al. (1983) Icarus 56:560 
koo- Fahr et al. (1991) J. Phys. Chem. 95:3218 
k^s- Tsang Sz Hampson (1986) J. Phys. Chem. Ref. Data 
15:1087 

A:q: Laufer et al. (1983) Icarus 56:560, reported in 
the NIST Chemical Kinetics Database V5.0. 

N.B. Laufer et al. presented tabular data; the NIST 
equation reproduces the Laufer et al. values to ~ 30% 
koo- Tsang (1989) Combust. Flame 78:71 

Zabarnick et al. (1986) J. Chem. Phys. 85:4373 

Tsang & Hampson (1986) J. Phys. Chem. Ref. Data 15:1087 

Becker et ah, (1989) Chem. Phys. Lett. 154:342 

Fahr et al. (1991) J. Phys. Chem. 95:3218 

Heinemann et al. (1988) Symp. Int. Combust. Proc. 21:865 

Monks et al. (1994) in preparation 






Notes: 


two body rates are in units of cm^ molecule"^ s“^ 

three body rates are in units of cm® molecule"^ s~^ 

if fc and /„ values are presented, then the factor F is given by 


logio F 


logic fc 



where M is the background atmosphere number density 



QUANTUM YIELDS IN METHANE PHOTOLYSIS 


The recent work of Mordaunt et al. [1993] by itself is incomplete, in that it does not provide sufficient information to 
unequivocally assign branching ratios (q) at Lyman a for all possible photolysis channels: 


CH4 + hi/ A 

b 

c 

d 


3CH2 + h + h 
1CH2 + H2 
CH + H2 + H 
CH3 + H 
ICH2 + H + H 


Furthermore, no measurements of ^CH2 or ^CH2 quantum yields are available and the data we currently have on CH3, CH, 
H2 , and H yields are contradictory. The experimental results on which we rely in our attempt to derive reasonable estimates 
of the branching ratios at 121.6 nm are: 

• net CH3 yield of 0.5 [Mordaunt et al. 1993]. 

• net H2 yield of 0.5 [Laufer & McNesby 1968]. Note: the experiment was actually performed at 123.6 nm with a reported 
quantum yield of 0.58, giving rise to the contradiction that qi, + qc = 0.58 but qd = 0.50. 

• net H yield of 1.0 [Slanger & Black 1982]. 

• net CH yield of 0.08 [Rebbert & Ausloos 1972]. Note: the measurements were at 123.6 nm (yield = 0.059) and 
104.8-106.7 nm (yield = 0.23); the estimate at Lyman a is from simple linear interpolation. 


These results, along with the requirement that the branching ratios for channels a-e sum to 1.0, are summarized by the 
following linear equation: 


1 1 

1 

1 

1 

qa 1.00 


0 0 

0 

1 

0 

qb 0.50 

CH3 yield 

0 1 

1 

0 

0 X 

qc = 0.50 

H2 yield 

2 0 

1 

1 

2 

0 

0 

H yield 

0 0 

1 

0 

0 

qe 0.08 

CH yield 

This matrix is singular, owing to the lack 

of constraints 

on *CH2 or ^CH 

2 yields, with channels a and e each yielding two 

hydrogen atoms. We choose to neglect channel e, which results in the overdetermined system 

1 1 

1 

1 




0 0 

0 

1 

9a 

CH3 yield 

(0.45 - 0.55) 

0 1 

1 

0 

X 96 

= H2 yield 

(0.45 - 0.55) 

2 0 

1 

1 

9c 

H yield 

(> 10) 

0 0 

1 

0 

9d 

CH yield 

(0.06 - 0.23) 


We solved this system by a least squares method, allowing for the variations in quantum yields given above, and searched 
for the “best” solution such that |1 — was minimized. It was necessary to renormalize the branching ratios to a sum 

of 1.0 once the “best” solution had been found; the loss of initial normalization stems from the dropping of one channel and 
from the “scatter” associated with the least squares solution technique. The result is 


9a 

= 0.212 

CH3 yield 

= 0.406 

96 

= 0.282 => 

H2 yield 

= 0.382 

9c 

= 0.100 

H yield 

= 0.930 

9d 

= 0.406 

CH yield 

= 0.100 


To extend these results to other wavelengths, we have adopted the following assumptions and constraints: 


• for A > 107 nm, the branching ratios for the various channels are constant with wavelength (the CH yield is 0.23). 

• for longer wavelengths, the CH quantum yield varies linearly with wavelength to reproduce the measured variation, 
while the CH3, H2 & H yields remain independent of wavelength (where energetically allowed). 

• 9 o(A > 132) = 0 and gc(A > 133) = 0. 

• photolysis branching ratios sum to 1.0. 

The laboratory measurements most needed are (1) the yield of ^CH2 or ^CH2 and (2) the yield of CH, both at Lyman a. 



DATA-MODEL COMPARISONS 


K profiles of the form 

with ^ = 0.5 (the o-subscript refers to a reference level) are used in generating the models 
shown in the top panels of the two adjacent figures, since this type of profile has been the 
one most commonly used in photochemical modeling studies of outer planetary atmospheres 
(typically blamed on Lindzen [1971] and Hunten [1975]). The fits are determined by the 
Ac < 140 nm channels, and it is apparent that the C 2 hydrocarbon abundances in these models 
are not consistent with the UVS data. In particular, the model lightcurves with Ac near 152 
nm exhibit more acetylene (C 2 H 2 ) opacity than is consistent with the data (these channels 
encompass the strong first vibrational member of the 3R-X Rydberg transition). At longer 
wavelengths, particularly in the egress case, the excess opacity in the models is from ethylene 
(C2H4). At shorter wavelengths, where ethane (C 2 H 6 ) is expected to contribute, the model 
lightcurves are in fairly good accord with the data (except for the egress case at Ac = 148.2 
nm, where the second member of the acetylene 3R-X system makes its appearance). Near 
the 0.5 mbax level, there is a large discrepancy with the IRIS-retrieved mixing ratios; this is 
attributed to K values in this vicinity too large to maintain the observed abundances. 

We subsequently defined AT-profiles characterized by a small number of layers over which K is 
held constant, and explored model fits by varying the width and placement of the constant-A" 
layers. The upper troposphere, lower stratosphere dynamics modeling of Conrath et al. [1991] 
suggests a sluggish overturning timescale on the order of 10® s, corresponding to an eddy 
mixing coefficient value of ~2 x 10^ cm^ s“^. Photochemical modeling of IRIS measurements 
of stratospheric acetylene and ethane emissions [Bezard et al. 1991; Romani et al. 1993] also 
strongly indicates a requirement for a sluggish lower stratosphere if the observed abundances 
are to be maintained. In view of this, the lowermost constant-A layer has been constrained to 
values on the order of 2 x 10^ cm^ s“^, although the extent of this layer has been varied. 

We found that a minimum of three constant- AT layers are required to give a good fit to the 
UVS data, with the maximum K values belonging to the middle layer centered near 10“^ mbar 
beneath the primary photochemical zone (where unit optical depth in methane at Lyman a 
occurs): 

I< = I<i, p> Pi 

= I<i X , Pl>p>P2 

= K 2 , P2> P>P3 

= K 2 X , P3>p>P4 

= Afg, P<P4 

where 71 = ln(Af 2 /Afi)/ ln(p 2 /pi) and 72 = ln(Af 3 //if 2 )/ In(p 4 /p 3 ). The placement of the rapid 
rise in K {pi > p > P 2 ) is chosen in part to reproduce a ratio of ethane-to- acetylene mixing 
ratios consistent with the Bezard et al. [1991] IRIS analysis and in part to avoid stifling the 
upward flow of CH4. Other types of profiles were explored {e.g., the profile form adopted by 
Yung et al. [1984] in their extensive Titan photochemical modeling study), but none were 
found to lead to fits as successful as models with K profiles defined by the above expression. 



The models shown in the lower panels of the adjacent figures illustrate the success (and model 
limitations) obtained with the latter form for K\ constituent opacities for these models are 
also shown. In the ingress case, models giving good fits for the Ac < 140 nm and Ac > 151 nm 
channels consistently failed to provide sufficient opacity at intermediate wavelengths, pointing 
to either an under- abundance of ethane or the presence of an unmodeled constituent. To 
gain an idea of the amount of C2H6 required to bring about agreement, the ethane mixing 
ratio profile was scaled by a constant factor; an increase in C2H6 abundances by a factor of 
4.0 in the model brings the model lightcurves into marginal alignment with the UVS data at 
Ac = 144.7 nm and Ac = 148.4 nm, but the agreement at Ac = 152.1 nm is now worsened. 
It is noteworthy, however, that the fit at shorter wavelengths is not compromised. It may 
be that our adoption of “local” conditions is inappropriate for the ingress data. If the rapid 
mixing in the middle stratosphere indicated by our models is a manifestation of advection in a 
meridional circulation system, then the photochemical stability of ethane would suggest that 
its abundance should be more reflective of global solar illumination conditions. 

The egress models are inherently close to such illumination conditions, and generating ethane 
abundances within the models sufficient to account for the opacities measured by the UVS is 
not difficult, using the A'-profile form with 3 constant layers. In the egress case, however, as 
in the ingress case, it is difficult to obtain consistent fits across the 140-153 nm interval, since 
ethane abundances yielding good fits at 140-148 nm result in too much total opacity at 152 
nm where the acetylene spectral signature lies. In the adjacent figure we show a compromise 
fit geared toward giving agreement with the IRIS-retrieved C2H6 mixing ratio at 0.7 mbar 
while maintaining a marginally acceptable fit to the UVS 152 nm lightcurve. A factor of ~3 
discrepancy with the IRIS acetylene results is noted. 

The K value in the middle layer is the minimum value needed in the models to prevent excessive 
C2H4 abundances vis-a-vis the data at wavelengths Ac > 153 nm. Why is an enhanced K 
necessary? - briefly, to transport C2H4 out of the primary photochemical zone where it is 
produced mainly by R3 (CH + CH4 — > C2H4 + H) to deeper levels where it is consumed in 
the three body reaction Rll (C2H4 -f- H + M — > C2H5 + M). The subsequent fate of C2H5 
is R12 (C2H5 + H — + CH3 + CH3), after which the methyl radical can combine with itself to 
form ethane (R6). Thus, enhanced K values in the vicinity of the primary photochemical zone 
(and beneath the nominal homopause) also contribute to ethane production. Even with K 
~ 10® cm^ s“^ in the middle layer, we have found it necessary to adopt a fa.st rate for Rll. 
Nevertheless, the inference of the existence of a region of enhanced eddy mixing at pressures 
beneath the nominal homopause is not critically dependent on the particular reaction rates we 
have adopted; it is based directly on the strong limits placed on ethylene abundances by the 
UVS data and on the displacement between C2H4 production and loss zones. 

At both ingress and egress, the placement of the “rapid-rise” in K is roughly consistent with 
the analysis of probable wave signatures in the RSS occultation data by Hinson and Magalhaes 
[1993]. However, we are not in a position to explain the large values of K in the middle layer 
of our models in terms of a particular mixing mechanism; in this work, K is simply a fitting 
parameter. However, it is our belief that, by and large, the chemistry is well accounted for, 
at least to the extent that our ignorance of dynamical conditions in Neptune’s stratosphere 
introduces greater uncertainties in modeling than remaining uncertainties in key reaction rates. 
If so, then our A-profiles yielding good agreement between models and data should provide 
useful constraints in dynamical modeling studies. 



To see if we could improve upon the model fits shown above, we have constructed models 
using piecewise constant K profiles partitioned by scale height above the tropopause. Given 
the rather large number of free parameters (more than can be uniquely specified by the data) , 
we have simply patterned these /{'-profiles on our previous results and attempted to refine the 
fits with the data. Typical results are shown in the adjacent figure. It heis not been possible (as 
yet) to improve upon the models based on K profiles of the 3-constant-layers variety, suggesting 
that more sophisticated models are not warranted by the data. The relative under- abundance 
of ethane is still apparent in the ingress model. The egress model was geared toward agreeing 
with both the IRIS C2H6 and C2H2 data points while remaining marginally consistent with 
the UVS data over the 140-153 nm range; the most obvious difficulty is with the UVS 151.9 
nm channel, owing to the abundance of acetylene required by the IRIS data. 

This analysis is continuing. 


Left Panels: thick solid curve, eddy mixing coefficient; solid curve, CH4 mixing ratios; long- 
dash curve, C2H6 mixing ratios; short-dash curve, C2H2 mixing ratios; dotted curve, C2H4 
mixing ratios. Scaled ethane mixing ratio profiles are shown in red. Altitudes referenced to 1 
bax radius of 24445 km (nominal atmosphere). 

Right Panels: solid curves, model Ijlo lightcurves; blue histograms, UVS data I /h ± <7. 
Lightcurves generated using scaled ethane mixing ratios are shown in red. Channel center 
wavelengths at 500 km minimum tangent ray height (MTRH) are indicated. 
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We have used a one-dimensional methane photochemical model 
to analyze Voyager observations of hydrocarbons and hazes in the 
stratosphere of Neptune. Voyager IRIS spectra provide informa- 
tkm about the global average C2H2 and C2H« mixing ratios for 
p > 0.1 mbar. The UVS lightcurves provide constraints on CH4 
and C2H4 in addition to C2H2 and C2H« but only at the solar 
occultation latitudes and for lower pressures. The model-predicted 
hydrocarbons arc very sensitive to the height profile of the eddy 
diffusion coefficient (/C). For both data sets K varying inversely 
with the atmospheric number density to some power produced 
poor results. Good agreement with the data requires that K be 
weak in the lower stratosphere (Af *• 2 x 1(P cm* see"* for p ^ 
2 mbar) but fairly vigorous in the upper stratosphere (If ^ S x 
10 * cm* sec~' for p% O.S mbar), i.e., a rapidly mixed upper 
stratosphere overlying a stagnant lower stratosphere with a rapid 
transition in between. The model C2H« and C2H2 mixing ratios 
are also sensitive to the reaction rate constants of C2H4 + H and 
CH] -E C2H3. Notably, we must use the present upper limit for 
the C2H4 -E H rate to best fit the model results to the observations. 
We are able to reproduce the IRIS C2H2 and C2H« emission features 
well, less so the UVS occultation lightcurves. Since the transport 
of C2H2, C2H«, and other hydrocarbons produced from methane 
photolysis out of the stratosphere is by ice haze formation and 
sedimentation, we compared model haze predictions to PPS and 
IRIS observations. For solar maximum fluxes (Voyager encounter 
conditions) the model mass production rate is 1 x 10“'* g cm* 
sec”'. C2H4 is the dominant haze component ( 75 %), with the 
remainder coming from C2H2 and C3 and C4 compounds. Balanc- 


ing the above haze production rate by the sedimentation rate for 
0.25-fim radius particles (upper limit to particle radius from PPS 
observations) yields a total haze column burden slightly above the 
PPS upper limit. However, lifetime analysis indicates that the 
model haze production rate should be averaged over solar mini- 
mum and maximum conditions. Under these conditions the model 
haze density is consistent with the PPS data. The predicted C4H2 
and C2H4 haze column densities are consistent with the lack of ice 
signatures in the IRIS spectra, o ms Academic Prm. inc. 


I. INTRODUCTION 

Methane (CH 4 ) and its photochemical products acety- 
lene (C 2 H 2 ) and ethane (C 2 HJ have been observed in 
emission in the infrared from Neptune’s stratosphere for 
over 10 years (see references in Orton and Appleby 1984, 
and Kostiuk el al. 1992). The retrieved mole fractions 
of these species are much greater than their respective 
saturation mixing ratios at the tropopause. This, coupled 
with observational evidence of a lower stratospheric haze 
(p > 5 mbar) on Neptune (see, e.g., Baines and Smith 
1990, Hammel el al. 1989), naturally led to the speculation 
that these and potentially more complex hydrocarbons 
were the source for the haze. 

An understanding of these hydrocarbons and hazes is 
important to an understanding of the carbon cycle in the 
stratosphere of Neptune. Carbon, in the form of methane. 
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